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Abstract 

The process of gluon-initiated double Higgs production is sensitive to non-linear interactions of the Higgs 
boson. In the context of the Standard Model, studies of this process focused on the extraction of the Higgs 
trilinear coupling. In a general parametrization of New Physics effects, however, an even more interesting 
interaction that can be tested through this channel is the tthh coupling. This interaction vanishes in the 
Standard Model and is a genuine signature of theories in which the Higgs boson emerges from a strongly- 
interacting sector. In this paper we perform a model-independent estimate of the LHC potential to detect 
anomalous Higgs couplings in gluon-fusion double Higgs production. We find that while the sensitivity to 
the trilinear is poor, the perspectives of measuring the new tthh coupling are rather promising. 



1 Introduction 



Measuring the couplings of the Higgs boson at the LHC is a difficult but important task. It will give 
crucial information to distinguish among different theoretical scenarios that can lead to a Higgs- 
like particle, and can thus shed light on the mechanism behind electroweak symmetry breaking 
(EWSB). In the Standard Model (SM), the request of perturbativity and unitarity up to Planckian 
scales fixes the strength of all the interactions of the Higgs boson in terms of its mass. Sizable 
modifications of the couplings can arise in weakly-coupled extensions, such as supersymmetry, 
through the mixing of the Higgs boson with new light states. In this case one expects to produce 
these new particles directly at the collider. A second compelling possibility is that the EWSB 
is triggered by new strong dynamics at the TeV scale, and a light Higgs emerges as the pseudo 
Nambu-Goldstone boson (pNGB) of a larger spontaneously broken symmetry [I]. In this case the 
modification of couplings is not necessarily accompanied by the presence of new light scalars and the 
direct manifestation of New Physics can be postponed to TeV energies. In addition to a modified 
pattern of linear couplings, this scenario predicts new non-linear interactions of the Higgs to the 
SM fields, which can lead to striking signatures at the collider and are a genuine feature of the 
underling strong dynamics. 

The problem of extracting the Higgs couplings by measuring its production and decay rates at 
the LHC has been studied at length in the literature, see for example Refs. [2]-[9]. The importance 
of a model-independent approach has been recently re-discussed, and a first estimate of the impact 



of the current LHC data on the Higgs parameter space has been performed in 10-19 



Aim of this paper is to study the effect of anomalous couplings in the process of gluon-initiated 
double Higgs production at the LHC, gg — > hh. Under the reasonable assumption of weak couplings 
to light fermions, this process proceeds through top quark loops and it is thus sensitive, in the first 
place, to the Higgs-top couplings. It also receives a contribution from the Higgs trilinear coupling, 



and for this reason has been studied in detail in the context of the SM 20-26 . In theories of 



New Physics, however, a much more interesting coupling that can be probed through this process 
is the non-linear interaction tthh. The latter is generally present in theories of composite Higgs, 
like for example the minimal models MCHM4 [27] and MCHM5 [28| , and gives a genuine signal of 
the Higgs strong interactions. As first noticed by the authors of Ref. [29], the presence of the new 
coupling can lead to a dramatic increase of the cross section. For example, enhancements larger 
than one order of magnitude are possible in the MCHM5, and even for (v/ f) 2 ~ 0.15, where / is 
the decay constant of the pNGB Higgs, the total cross section doubles compared to its SM value. 
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Given that the deviations due to Higgs compositeness are usually much milder, gg — > hh seems an 
extremely favored channel which is worth investigating. 

In this paper we derive a first quantitative assessment on the detectability of the anomalous 
coupling tthh in the process gg —> hh. In section 2 we briefly summarize the parametrization of 
the couplings of a generic Higgs- like scalar which we adopt. In section 3 we analyze the gg — > hh 
process in the presence of modified Higgs interactions. In particular we study the dependence 
of the cross section on the various couplings and show that there is high sensitivity to the new 
tthh coupling, much larger than that on the trilinear self-interaction. For our analysis we wrote 
a dedicated computer code which computes the exact 1-loop matrix element for single and double 
Higgs production via gluon fusion as a function of the relevant couplings. The code has been 
implemented as one of the available processes of the event generator ALPGEN 1 30 1 and will be 
made public with its next official release. We then discuss two of the most promising decay channels 
of the Higgs pair: hh — > WWjj — > Ivjj^^ in the case in which the Higgs has a suppressed single 
coupling to the top (fermiophobic limit), and hh — > bbj'j in the case in which the linear couplings 
are SM like. For the latter case, we follow the strategy proposed in Ref. [26] and in section 4 we 
perform a first collider study to estimate the exclusion and discovery limits on the anomalous tthh 
coupling. To compute the SM background cross section we use the results of [26] with updated b 
and 7 efficiencies and rejection factors. We collect the results of our collider study in section 4.1. 
Finally, conclusions are reported in section 5. 

2 General parametrization of the Higgs couplings 

In this section we introduce the general parametrization of Higgs couplings that will be used in this 
paper. 

The most general effective Lagrangian that parametrizes the interactions of a Higgs-like scalar 
at low energy has been discussed in 31 and extended in [IT]. Under the assumption of custodial 



symmetry, the Nambu-Goldstone bosons associated to the electroweak symmetry breaking can be 
described as the coordinates of the coset SU{2)i x SU(2)r/ SU{2)y ~ 50(4)/ SO(3). They can be 
conveniently parametrized by the 2x2 matrix 

Z = exp{ia aX a (x)/v) , a = 1,2, 3 (1) 

where a a are the Pauli matrices and v = 246 GeV. At energy scales much below possible new 
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physics states, the effective Lagrangian describing a light Higgs h has the form 

C = ^Tr(D^D^) (l + 2a £ + . . .) + \{d,hf - - d 3 \ (^f\ h 3 + . . . 

(2) 

- m t q^EntR ^1 + c t ~ + c 2 ^ + • • ■) ~ m b Ql^or f 1 + c b ~ + . . . J + h.c. , 

where qi = (£z,j &l) and a, ct t b, c 2 and c?3 are the numerical coefficients that parametrize the Higgs 
couplings. The dots stand for terms which are not relevant for double Higgs production via gluon- 
fusion. In particular, we assume that the strength of single interactions of the Higgs to the fermions 
is not extremely enhanced compared to its SM value, so that the contribution of the light fermions 
and the bottom quark to double Higgs production can be safely neglected. The coupling to the 
bottom, Cf,, is relevant only in the Higgs decay, and we will set 

c t = c b = c (3) 

for simplicity in the following. We also neglect ggh and gghh local interactions which can be 
generated by new heavy states at 1-loop level, as for example scalar or fermionic partners of the 
top quark. 

In this analysis we will freely vary the parameters that appear in the effective Lagrangian. In 
specific models, however, they can be related to each other. For example, the SM Lagrangian is 
obtained for 

a = c = d 3 = 1 , c 2 = . (4) 

A class of theories that we will consider in the following are the composite Higgs models based on 



the symmetry pattern SO(5)/SO(4) (27 28 . In these models the parameter a is given by 

«=\/W, (5) 

where £ = v 2 /f 2 and / is the Nambu-Goldstone decay constant. The values of c, c 2 and d 3 depend 
on which SO (5) representation the fermions are embedded in. For the two minimal choices of 
fermions in the spinorial (MCHM4 |27|) and fundamental (MCHM5 [28]) representations one gets 



c = d3 = \/l — C2 = — - , MCHM4, spinorial representation , (6) 

1 - 2£ 

c = d 3 = , C2 = — 2£ , MCHM5, fundamental representation . (7) 

yl - £ 

Equations ^ , Q and ([7]) account for the value of the Higgs couplings as due to the non-linearities 
of the chiral Lagrangian. The exchange of new heavy particles can however give further corrections 
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Figure 1: Feynman diagrams for double Higgs production via gluon fusion (an additional contribution comes 
from the crossing of the box diagram). The last diagram contains the new non- linear Higgs interaction tihh. 



to these expressions. In the following we will neglect these effects since they are parametrically 
subleading [32], although they can be numerically important when the top or bottom degree of 
compositeness becomes large 33 . This is especially justified considering that in minimal composite 
Higgs models with partial compositeness these additional corrections to the couplings do not affect 
the gg —> h rate because they are exactly canceled by the contribution from loops of heavy fermions, 



as first observed in Refs. [34, 32 and explained in Ref. 33 . For double Higgs production we 
expect this cancellation to occur only in the limit of vanishing momentum of the Higgs external 
lines. In general, numerically important contributions might come from light top partners (light 
custodians). In models with partial compositeness, where the dominant contribution to the Higgs 
potential comes from top loops, the presence of light fermionic resonances is essential to obtain 
a light Higgs [28 35 . In particular, mh ~ 120 — 130 GeV requires top partners around or below 
1 TeV. It would be interesting to analyze in detail their effects on double Higgs production. 



3 Double Higgs production via gluon fusion 

In the scenario we are considering, the leading-order contributions to the process gg —> hh come 
from Feynman diagrams containing a top-quark loop. The three relevant diagrams are shown 
in Fig. [TJ and can be computed by using the results of Ref. [21 J . We have implemented the 
automatic computation of the matrix element as one of the processes of the ALPGEN MonteCarlo 
generator |30| . The code will be made public with the next official release of ALPGEN, and it allows 
one to compute the total cross section and differential distributions, as well as to generate events for 
an arbitrary choice of the Higgs couplings c, d^, ci- The validation of the code has been performed 
by means of an independent C++ program linked to the QCDLoop [36] and to the LHAPDF 
routines |37|. All the results reported in the following have been derived by use of the ALPGEN 
matrix element calculation with CTEQ61 parton distribution functions and renormalization and 
factorization scales Q = m{hh). The top quark mass has been set to mt = 173 GeV. 

The amplitude of each diagram in Fig. [T] is characterized by a different energy scaling at large 
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invariant masses 



m(hh) 3> mt,mh. One has 



a 2 m t 

An ~ c a s — 



Aa ~ c c?3 a. 



mr m. 



v 2 s 



log 



mr 



Anl ~ c 2 a s 



log 



mr 



+ Z7T 



+ ITT 



(8) 
(9) 

(10) 



where An , A a ar e the amplitudes of respectively the box and the triangle diagram with the Higgs 
exchange (first two diagrams of Fig. [I]), while A Anl denotes the amplitude of the diagram with the 
new non-linear interaction tihh (last diagram of Fig. [TJ. At large s the box and the diagram with 
the new vertex dominate, while the triangle with Higgs exchange gives its largest contribution near 
threshold. For SM values of the couplings there is a destructive interference between Aa and An, 
so that decreasing the trilinear coupling d 3 leads to a softer distribution, while increasing it makes 
the suppression of the cross section near threshold even stronger. On the other hand, since AahI 
and An have similar energy scalings (the log enhancement of AauI becomes important only at very 
large Vs where the gluon pdfs are small), the main effect of their interference is on the total cross 
section, with little modification of the m(hh) distribution. These behaviours are clearly visible in 
the distributions shown in Fig. [2] for m^ = 120 GeV at 14 TeV. 

The cross section depends on the couplings as a quadratic polynomial in the variables c 2 , c 2 and 
cd 3 , associated respectively to the three diagrams of Fig. [TJ It can thus be conveniently expressed 
by the formula 



° (PP ~~ ^ hh) 



4 + («c 2 ) 2 + (f3cd 3 ) 2 + A lC2 (ac 2 ) + A 2 (ac 2 ) (f3cd 3 ) + A 3 c 2 (f3cd 3 ) 



fill 



where the value of the (real) coefficients a, a, (3, A\, A 2 , A 3 has been extracted by fitting the results 
of a Montecarlo integration, and is reported in Table [l] (at LO in a s ) for m^ = 120, 125 GeV at 8 
and 14 TeV. In the above parametrization, the coefficients a and f3 measure the sensitivity of the 
cross section on the parameters c 2 and (cd 3 ), relative to c 2 . One can see from Table [l] that the 
dependence on c 2 is significant while the one on (cd 3 ) is rather mild. This can be tracked back to 
the additional factor (m^/s) in the amplitude of the triangle diagram which carries the dependence 
on (cd 3 ), see eq.Q. This factor leads to a suppression at large s and thus, because of the kinematic 
threshold s > 4m 2 , to a reduction of the sensitivity on (cd 3 ) of the total cross section. 

We see in Table[T]that increasing the LHC center-of-mass energy from 8 TeV to 14 TeV increases 
the LO total cross section (hence the coefficient a) by a factor ~ 5.2, while the relative strength 
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Figure 2: Invariant mass distribution of the two Higgs bosons in pp — > hh at the LHC (14 TeV) for 
to/j = 120 GeV. The various curves correspond to different choices of Higgs couplings: c = d,3 = 1, C2 = 
(SM couplings, solid blue curve), c = d 3 = 1, c 2 = —1 (dashed purple curve), c = 1, d 3 = c 2 = (dotted 
yellow curve). The dot-dashed green curve shows the distribution obtained in the approximation of infinite 
top mass with SM couplings. All curves have been normalized to unit area. The corresponding LO total 
cross sections are 15.2 fb (solid blue curve), 253 fb (dashed purple curve), 31.6fb (dotted yellow curve). 
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Ax 


A 2 


A 3 


120 GeV 


14 TeV 


151.3 fb 


0.453 


0.164 


-1.86 


-1.77 


1.66 


8 TeV 


32.6fb 


0.474 


0.178 


-1.89 


-1.78 


1.68 


125 GeV 


14 TeV 


144.6 fb 


0.457 


0.169 


-1.85 


-1.79 


1.68 


8 TeV 


30.5 fb 


0.475 


0.185 


-1.89 


-1.79 


1.70 



Table 1: Coefficients for the fit of eq.( 11 1 of the total LO pp —> hh cross section via gluon fusion at the LHC. 



among a, f3, A\,A2, A3 varies by less than ~ 15%. When the Higgs mass is varied from 120 GeV to 
125 GeV, a decreases by ~ 5 — 7%, while the other coefficients change by less than ~ 1 — 3%. The 
left plot of Fig. [3] illustrates how the total cross section changes when varying individually (I3 and 
C2, while fixing the other couplings to their SM value. In the vicinity of the SM point, decreasing 
(increasing) or c<i leads to an enhancement (reduction) of the total cross section, with a much 
stronger dependence on C2 than on d^. The right plot in the same Figure shows how the cross 
section varies with the Higgs mass for the SM choice of couplings. The solid curve corresponds to 
the full one-loop matrix element calculation, while the dot-dashed curve is obtained by taking the 
limit of infinite top mass. As previously noticed [23], Q this approximation is reasonably accurate in 



1 The infinite top mass limit is also discussed by the authors of Ref. 20 , although the contribution of the triangle 
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Figure 3: Left plot: total cross section in SM units as a function of C2 = A (solid blue curve) and ds = 1 + A 
(dashed purple curve), with the other Higgs couplings set to their SM values and = 120 GcV. Right plot: 
LO total cross section (in fb) in the SM as a function of rrih as computed by means of the full one-loop 
matrix element (solid blue curve) and the infinite top mass approximation (dot-dashed green curve). In both 
plots the LHC center-of-mass energy has been set to 14 TeV. 



the case of the total cross section, but completely fails to reproduce the correct m(hh) distribution, 
as illustrated by the corresponding curve in Fig. [2] 

We have seen that modified Higgs couplings, in particular a non- vanishing tthh interaction, can 
lead to a strong enhancement of the total pp — > hh cross section. However, in order to determine 
the signal yield at the LHC in a given final state one has to take into account also the change in the 
Higgs decay branching ratios. In our case, these latter depend only on the ratio of the parameters 
c and a. We will consider two illustrative situations: i) the case in which the branching ratios are 
similar to the SM ones; ii) the case in which the couplings of one Higgs boson to two fermions 
are suppressed (fermiophobic limit). The first situation is realized in models where all single Higgs 
couplings are rescaled by the same factor (as in the MCHM4, see eqs.Q, ([6])), or their shift from 
the SM value is small (as in composite Higgs models with small £). In this case the studies of 



Refs. (23-26 suggest that the most favorable final state at the LHC for a light Higgs boson is 
hh — > 7766. On the other hand, if single Higgs couplings to fermions are suppressed (while the 
strength of tthh can still be sizable), the dominant decay mode for a light Higgs becomes h — > WW. 
The channel hh — > WWWW has in this case the largest rate and should be visible in final states 
with two or three leptons. The 77 branching ratio is also strongly enhanced, so that in this case 
hh — > WW'jj also seems a promising final state. In particular, the request of one lepton from the 



diagram seems to have been accidentally omitted in their Fig. 2. 
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Figure 4: Value of the branching ratio BR(hh — > 7766) (on the left) and BR(hh — > WWjj) (on the right) 
in SM units as a function of the ratio of Higgs couplings c/a. The dots show the prediction in the MCHM5, 
where c/a = (1 — 2£)/(l — £), for various values of £. In both plots the Higgs mass is set to rrih = 120 GeV. 



decay of the W pair should be sufficient to reduce the background and lead to a clean signature at 
the LHC. 

Figure [4] illustrates how the branching ratios BR(hh — > 7766) and BR(hh — > WW^y^f) vary 
with c/a. Strong enhancements compared to the SM prediction are possible for BR(hh — > WW-fj) 
in the fermiophobic limit c — > 0. A fermiophobic composite Higgs can for example arise in the 
MCHM5 for £ — > 1/2, see eq. Q. Although the point (a = 1, c = 0) has been excluded at 
95% CL in the range = 110 — 192 GeV by the combination of all CMS searches 38 , the one 



predicted by the MCHM5 for £ = 1/2 (a = 1/V2,c = 0) is still allowed for m h ~ 125 GeV and 
in fact could better explain the pattern of observed enhancements in the various 77 categories of 
the CMS analysis 39 . Figure [5] shows the final yield per fb _1 predicted in the MCHM5 in the 
two final states hh — > 77^6 and hh — > WWjj — > Ivqq^ as a function of £ for = 120 GeV. 
The rate has been computed using the cross section for gg — > hh at LO in a s {i.e. no iT-factor 
is included) given by eq.(ll) and Table [T] As expected, for small values of £ the most promising 
channel is 7766, whose rate can be significantly enhanced compared to the SM expectation. At 14 
TeV, for example, even for £ = 0.1 the signal yield more than doubles. This large sensitivity to 
small values of £ = (v/ f) 2 shows that double Higgs production via gluon fusion is an extremely 
powerful process to probe the Higgs compositeness at the LHC. Still, the difficulty of isolating the 
77&5 signal from the background will require large integrated luminosities and will be possible only 
in the high-energy phase of the LHC. In the fortunate situation in which the Higgs is fermiophobic, 
on the other hand, the enhancement of the hh —> WW'yy — > Ivqq^ final state is so large in the 
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Figure 5: Signal yield per fb _1 predicted in the MCHM5 for the two final states hh — > 77&B (solid curves) 
and hh — > WWjj — > Ivqepj^j (dashed curves) as a function of £. The thick (thin) curves correspond to the 
LHC with 14 TeV (8 TeV) center-of-mass energy. The Higgs mass is set to m/j = 120 GeV. The rate has 
been computed using the cross section for gg — > hh at LO in a s (i.e. no K- factor is included) given by 
eq.([TT]) and Table [l] 



MCHM5 that a first preliminary observation of the signal might be possible at the 8 TeV LHC. 
At this energy, for £ = 1/2, m h = 120 GeV and 20 fb _1 of integrated luminosity the MCHM5 
predicts ~ 15 signal events before cuts. Considering that the SM background is expected to be 
rather small for a final state with two photons and one isolated lepton, this number of events might 
be sufficient to establish the observation of the signal. For the same value of c.o.m. energy and 
integrated luminosity, the MCHM5 predicts ~ 42 and ~ 27 signal events (before cuts) respectively 
in hh — > 4W — > X^X^vvAq (two same-sign leptons) and hh — > 4W — > ?,l2>uqq. These high rates 
suggest that it might be possible to distinguish the hh — > AW signal over the SM background even 
at 8 TeV. 

The results discussed in this section are rather encouraging, and show that double Higgs produc- 
tion can be an important process to extract or constrain the tthh interaction and, more in general, 
to probe the Higgs compositeness. However, a more robust assessment of the LHC sensitivity in this 
sense requires a dedicated analysis of each of the relevant final states and a careful estimate of the 
background. In the next section we will focus on the 7766 channel and use the studies of Ref. 
to get a first determination of the precision which can be obtained on C2 at the 14 TeV LHC. 



26 



We included a Jf-factor = 2 in the estimate, which is the value obtained in Ref. 22 at 14 TeV, assuming that 



similar result also applies at 8 TeV. For mh = 125 GeV the number of signal event is ~ 10. 
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4 Analysis of the 6677 channel 



The analysis of the gg — >• hh — > 6677 process performed in Ref. 26 aimed at measuring the Higgs 
trilinear coupling, and assumed SM values for the other couplings. In this section we make use of 
the results of [26] to estimate the LHC sensitivity on the tihh non-linear interaction. 

We assume a center-of-mass energy y/s = 14 TeV and set = 120 GeV. Given the value of the 
signal rate and the large SM background, an analysis of 6677 at 8 TeV seems rather challenging and 
for this reason it will not be considered here. We expect our results at 14 TeV to be representative 
of what the LHC sensitivity will be in its future high-energy phase, even if the actual value of 
c.o.m. energy turns out to be different (on the prospects of 14 and 13 TeV see for example |40|). 
The analysis is performed at the parton level: signal events are generated by means of our im- 
plementation of double Higgs production via gluon fusion in ALPGEN, while the computation of 



the background processes is taken from 26 . We take the NLO QCD corrections to gg — > hh into 
account by multiplying the LO cross section by a factor K = 2 (22] , [^] while we neglect the smaller 
contribution of the vector-boson-fusion process to double Higgs production. To ensure an effective 
suppression of the background, we select events with two photons and two 6-jets. Two photon tags 
are necessary to suppress the huge QCD background. On the other hand, the optimal number of 



6-tags depends on the strength of the signal and on the collider energy. In Ref. 26 only one 6-tag 
was required for the analysis at the 14 TeV LHC. This choice was motivated by the necessity of 
preserving the small SM signal as much as possible, at the expense of having to cope with a larger 
background. In our case the signal cross section is much larger than the SM one in a sizable part 
of the parameter space. Furthermore, recent studies have shown that a larger 6-tagging efficiency 
than that assumed in Ref. |26| is possible at the LHC [41] , maintaining an acceptable rejection rate 
on jets. In particular, we adopt the following conservative estimates for the efficiencies (e) and fake 
rates (r) for 6-tagging 41 and photon reconstruction |42| : [^| 



e 7 = 0.8, r 7 = 2500, e 6 = 0.7, r c ^ b = 5 , r^ b = 25 . (12) 

This allows us to require two 6-tags, thus obtaining a stronger suppression of the background at 
the price of an affordable reduction of the signal. We find that including the additional category 
of events with one 6-tag in the analysis does not sensibly improve our results. 

3 Notice that the authors of Ref. [26| used a factor K — 1.65 for the LHC at 14 TeV, which is the one appropriate 
for their choice of the renormalization and factorization scale Q — mu- As previously discussed, we set instead 
Q — m(hh). 

4 The authors of Ref. |26| use instead tb = 0.5 and fake rates r c -»b = 13 and r^b = 140. 
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used cuts 


6677 


CC77 


bbjj 


cc-yj 


jjn 


bbjj 


ccjj 


1333 


3333 


eq. (13) 
+ eq. (14) 


0.056 
0.0060 


0.42 
0.0215 


65 
8.28 


250 
17.0 


11 

0.84 


2.5 xlO 4 
4520 


2.5xl0 4 
4520 


7700 
364 


5xl0 6 
4xl0 5 


+ tags 


0.0019 


5xl0" 4 


0.0013 


2xl0" 4 


9xl0" 4 


4xl0~ 4 


3xl0~ 5 


2xl0" 4 


lxlO" 4 



Table 2: Cross sections (in fb) of the main QCD backgrounds to hh — > bbjj at the 14 TeV LHC. The values 
of the cross sections after the cuts of eqs. ( 13 ) and ( 14 ) are taken from (26 , and include an additional rescaling 



factor 1.3 introduced to take into account the increase that can possibly come from NLO corrections. The 
last line reports the value of the cross sections after the inclusion of the 6-jet and photon tagging efficiencies 



of eq. (121. 



Following 26 , we impose the set of kinematic cuts 

p T {b)> 45 GeV, |t/(6)|<2.5, AR(b, b) > 0.4 , 
m h - 20 GeV < m(bb) <m h + 20 GeV , 

p r (7)>20GeV, |t/( 7 )|<2.5, Ai?( 7 , 7) > 0.4 , (13) 
m h - 2.3 GeV < 771(77) < ™h + 2.3 GeV , 
Ai?( 7 ,6) > 0.4, 

which have a high selection efficiency on the signal and ensure that the 66 and 77 invariant masses 
are reconstructed in the given windows around the Higgs mass. To take into account the detector 
resolution, we have assumed a 79% efficiency for the reconstruction of the bb pair and a 79% 
efficiency for the reconstruction of the 77 pair in the signal, as done in Ref. |26|. After the above 
cuts, the most important irreducible backgrounds come from the 6677 continuum, and potentially by 
single Higgs production in association with two 6-quarks (h(—> 77)66) or two photons (h(— > 66)77). 
The reducible backgrounds are QCD processes (cc 77 , 6677, cc 7 j, 66jj, ccjj, njj, ijjj, jjjj) or 
single-Higgs production processes (hjj and /177) where some of the jets fake one or more 6-quarks 
or photons. While the QCD backgrounds are universal and as such are not affected by New Physics, 
those coming from single Higgs production depend on the value of the modified Higgs couplings. 
In the SM, all single Higgs processes are much smaller than the double Higgs signal after the cuts 



of eq.(13), and can be safely neglected |26|. This approximation is still valid in our context, and 



for this reason we will not include these backgrounds in our analysis. The list of relevant processes 



and their cross sections after the cuts of eq. (13) (without including 6 and photon reconstruction 
efficiencies) is reported in Table [2] 
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A further suppression of the background can be obtained by exploiting the particular topology 
of the signal, where the two Higgs bosons are produced back to back in the center-of-mass frame. 
We thus select events where the 77 pair has a small opening angle, while the minimal angular 
separation between a 6-jet and a photon is large. We require [26] : 

Ai?( 7 , b) > 1.0 , A,R(7, 7) < 2.0 . (14) 

In most of the parameter space these cuts imply a moderate reduction of the signal (20% — 40%), 
while the total background is suppressed by one order of magnitude. The corresponding background 
cross sections are reported in the second line of Table [2} 

Finally, the reducible backgrounds are drastically suppressed once the efficiencies for recon- 



structing two photons and two 6-jets of eq.(12) are included. The resulting final cross sections are 



shown in the last line of Table [2] The corresponding total background cross section is = 5.5 ab. 



After performing the kinematic cuts of eqs. (13), (14) and including the efficiencies for the 
reconstruction of the 77 and bb pairs, the tagging efficiencies (12), and the -RT-factor, the signal rate 
(r s = cr(pp — > hh) x BR(hh — > 77W)) at = 120 GeV is well approximated by the formula 

BR(hh -> 7766) 



BR{hh — > ^bb)sM 



(49.3ab) c| + (0.407c 2 ) + (0.101 cd 3 f - 1.76 c 2 (0.407 c 2 ) 

(15) 

1.82 (0.407 c 2 ) (0.101 cd 3 ) + 1.72 c 2 (0.101 cd 3 ) 



For the SM case (c = d 3 = 1, c 2 = 0) we find r s = 4.9 ab. Notice that, compared to the fit of 
eq.( 11 ) at 14 TeV (see Table [T]), the cuts have further weakened the dependence of the cross section 



on d 3 , since their efficiency is smaller for events with low m(hh) invariant mass. 
4.1 Results 

Using the signal and background rates derived above, we can estimate the sensitivity of the 
14 TeV LHC on pp — > hh —> 6677 for = 120 GeV. Since large luminosities are typically needed 
to distinguish the signal from SM background, we expect that by the time the analysis of double 
Higgs production is performed, the Higgs branching ratios to 77 and bb and the linear couplings 
o, c are known with good accuracy. Double Higgs production can thus be used to extract (or set 
limits on) the couplings c 2 and c^. 

Figure [6] shows the luminosity required to discover the signal as a function of c, c 2 and ^3, 
assuming that the branching ratio BR(M — > 7766) has the value predicted in the SM.^The plots 



The definition of discovery luminosity and the details of our statistical analysis are discussed in the Appendix. 
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Figure 6: Isocurves of discovery luminosity (in fb" 1 ) at the 14 TeV LHC in the plane (02,^3) for c = 1 
(on the left) and in the plane (c, C2) for = c (on the right). Outside each contour, the pp — > hh — > 7766 
signal can be discovered with the corresponding integrated luminosity. In both plots the Higgs mass is set 
to rrih = 120 GeV and the Higgs decay branching ratios are fixed to their SM values. See the Appendix for 
the definition of discovery luminosity. 



on the left and on the right show the luminosity contours respectively in the plane (02,^3) for 
c = 1, and in the plane (c, C2) for d% = c. As expected, the sensitivity on c and C2 is stronger than 
that on the Higgs trilinear coupling ^3. In particular, while a discovery in the SM would require 
at least 1200 fb~ , we find that much lower luminosities are sufficient even for moderately small 
values of C2. Figure [7] shows the corresponding discovery luminosity in the composite Higgs models 
MCHM4 and MCHM5 as a function of £. We find that values of £ as small as 0.15 can be probed 
with 300 fb -1 of integrated luminosity. Compared to other processes like double Higgs production 
via vector boson fusion [31] , these results show that gg — > hh can be extremely powerful to study 
the non-linear couplings of a composite Higgs and thus probe its strong interactions. 

Once a discovery is established, one can measure the couplings C2 and d% by using the value of 
c and of the Higgs branching ratios determined in single-Higgs processes. The left plot of Fig. [| 
shows the region of 68% probability in the plane (C2, d^) with 300, 600 and 1200 fb -1 (light, medium 
and dark blue regions) obtained by injecting the SM signal (c = 0^3 = 1, C2 = 0). H In this case 
the precision on C2 is poor even with 1200 fb^ 1 , while d% is basically unconstrained. A much 
6 That is: the rate of observed events is assumed to be that predicted in the SM with = 120 GeV. 
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Figure 7: Discovery luminosity at the 14 TeV LHC in the MCHM4 (dashed blue curve) and MCHM5 
(continuous purple curve) as a function of £. The Higgs mass is set to = 120 GeV. See the Appendix for 
the definition of discovery luminosity. 
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Q c 2 

Figure 8: Regions of 68% probability in the plane (02,^3) obtained with 300 (light blue area), 600 (medium 
dark blue area) and 1200 fb _1 (darker blue area) of integrated luminosity. On the left: injected signal is the 
SM (c = d 3 = 1, c 2 = 0); On the right: injected signal is the MCHM5 with £ = 0.3 (c = d 3 = 0.48, c 2 = -0.6). 
Both plots are obtained by fixing the Higgs branching ratios and the coupling c to their corresponding 
theoretical values, assuming that they can be measured accurately from single-Higgs processes. 
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Figure 9: Precision on the parameter £ in the MCHM5 which can be obtained at the 14 TeV LHC for 
mh — 120 GeV from the analysis of pp —> hh —> 7766. For each injected value £ t h, the solid, dashed and 
dot-dashed curves (red, orange and yellow regions) show the 68% probability interval which is expected on 
the measured value £, exp respectively with 600, 300 and 100 fb _1 . The gray regions delimited by the solid, 
dashed and dot-dashed vertical lines indicate the values of £, respectively with 600, 300 and 100 fb _1 , for 
which the signal rate is too small to make a discovery, (see Fig. In particular, for £ — > 0.5 the h — > bb 
decay rate vanishes (the Higgs becomes fermiophobic) , and the signal cannot be distinguished from the SM 
background. The curves are derived by assuming that the 77 and bb branching ratios are determined with 
good accuracy from single- Higgs production. 

more precise determination of C2 can be obtained if its value is non- vanishing. The right plot of 
Fig. [8] shows the case in which the injected signal is the MCHM5 with £ = 0.3, corresponding to 
(c = = 0.48, C2 = —0.6). It assumes that the Higgs branching ratios and the coupling c = 0.48 
predicted by this model have been measured with good accuracy in single-Higgs processes. We find 
that with 300 fb _1 the coupling ci can be determined, up to a discrete ambiguity, with a precision 
of ~ 15 — 20%. On the other hand, even in this case remains largely unconstrained with our 
analysis. Finally, Fig. [9] shows how precisely the parameter £ can be determined in the MCHM5 
through gg — > hh — > 7766 by making use only of the value of the decay branching ratios determined 
in single-Higgs processes (that is: without fixing c to its measured value in the fit). For each 
injected value ^th, the solid, dashed and dot-dashed curves (red, orange and yellow regions) show 
the 68% probability interval which is expected on the measured value £ exp respectively with 600, 
300 and 100 fb _1 . For example, with 300 fb _1 of integrated luminosity, £ = 0.2 can be measured 
with a precision of ~ 30%, which is in agreement with Fig. [8] and shows that the precision on £ is 
driven mainly by that on C2. The solid, dashed and dot-dashed vertical lines indicate the range of 
£, respectively for 600, 300 and 100 fb _1 , for which the expected signal yield is sufficiently large to 
establish a discovery (see Fig. [7]). 
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5 Conclusions 



The discovery of a light Higgs-like scalar at the LHC will mark a first important step forward in 
our comprehension of the electroweak symmetry breaking mechanism. Precise knowledge of the 
strength of its interactions with SM fields can shed light on the origin of the Higgs boson and 
indicate if the new dynamics at the electroweak scale is weakly or strongly interacting. New strong 
dynamics can form the Higgs as a bound state and solve naturally the hierarchy problem of the 
Standard Model. In this case the Higgs boson itself interacts strongly at large energies due to its 
modified linear couplings to SM fields and the existence of new non-linear interactions. 

In this paper we have performed a first model-independent study of double Higgs production 
via gluon fusion, gg — > hh, and we have shown that its cross section is greatly enhanced by the non- 
linear interaction tthh. Such new vertex gives a contribution to the scattering amplitude of tt — > hh 
that grows with the energy, A(tt — > hh) ~ (Emt)/v 2 , and as such it is a genuine signature of the 
underlying strong dynamics. In the process gg — > hh the tthh vertex mediates a new diagram 
containing a top-quark loop which grows logarithmically at high energies and does not lead to 



a violation of perturbative unitarity (see eq.(10)). However, it does lead to a strong numerical 
enhancement of the cross section compared to the SM, as first noticed by the authors of Ref. |29| in 
the context of the MCHM composite Higgs models. The origin of this enhancement can in part be 
traced back to the sizable destructive interference which occurs in the Standard Model between the 
box and the triangle diagram with Higgs exchange. A similar cancellation has been found to take 
place at large s in gg — > ZZ [43] and gg — > WW [44], and interpreted as a relic of the cancellation 
dictated by unitarity between the energy-growing amplitudes in the sub-process tt —> VV . On the 
other hand, in the case of double Higgs production the cancellation is in the total cross section (i.e. 
not necessarily at large s). Furthermore, none of the diagrams which contribute to tt — > hh in the 
SM grows with the energy, so that the cancellation between the box and the triangle diagram in 
gg —> hh should be rather seen as a numerical accident, not driven by the unitarity of its subprocess. 

The strong enhancement of the gg — > hh cross section makes this process quite powerful to 
measure or constrain the strength of the tthh interaction, which we have denoted as C2- This 
should be compared with the much weaker sensitivity on the Higgs trilinear interaction, d%, which 
makes the extraction of this latter coupling extremely challenging at the LHC. In particular, the 
weaker dependence on d^ follows from an extra suppressing factor (m^/s) carried by the triangle 
diagram with Higgs exchange, which thus contributes mainly at threshold. In order to assess 



the LHC precision on C2 we have made use of the results of Refs. 23-26 , where double Higgs 
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production via gluon fusion was studied in the context of the SM and several Higgs decay channels 
were investigated. In a generic scenario of New Physics, what is the best final state largely depends 
on the value of the Higgs decay branching ratios. In particular, enhanced branching ratios can 
combine with the increase in the double Higgs production cross section and lead to dramatic effects 
at the LHC. For example, in the (fermiophobic) limit in which the linear couplings of the Higgs 
to the SM fermions are suppressed (c — > 0), the h — > WW and h — > 77 branching ratios can be 
sensibly enhanced compared to their SM values. In this case the final states hh — > WWWW and 
hh — > WW77 seem to be extremely promising, and might be visible even at the 8 TeV LHC. At 
this energy and with L = 20 fb , for example, we find that for rrih = 120 GeV the MCHM5 at 
£ = 0.5 predicts ~ 15 signal events in hh — > WWjj — > Ivqq^ (see Fig. [5]). 

If the Higgs decay branching ratios do not differ much from their SM values and the Higgs 
boson is light, the most powerful final state should be hh — > 6677, as suggested by the analysis 
of Ref. [26]. Even in this case, the signal rate can be significantly enhanced compared to the SM 
prediction if C2 is not too small. At 14 TeV, for example, the signal rate predicted in the MCHM5 
with £ = 0.1 is larger than the SM one by more than a factor two. This sensitivity on small values 
of £ = {v/f) 2 shows that double Higgs production via gluon fusion is an extremely powerful process 
to probe the Higgs compositeness at the LHC. In order to estimate the LHC sensitivity on C2 in 



a model-independent way, we have followed the strategy proposed in Ref. 26 and performed a 
Montecarlo study of pp — > hh — >• 6677. We have computed the SM background cross section by 
using the results reported in I26J and rescaling them to take account of updated b and 7 efficiencies 
and rejection factors. The results that we obtained are quite encouraging. With L = 300 fb -1 the 
14 TeV LHC can probe values C2 < —0.2 and C2 > 0.8 if c, d% ~ 1 (see Fig. [6]). In the case of the 
MCHM5, an integrated luminosity of 300 fb -1 is sufficient to discover a signal with £ > 0.15 (see 
Fig. [7]). In general, once the signal can be statistically distinguished from the background and a 
discovery is made, the value of C2 can be extracted with good accuracy. For example, we find that 
by injecting a signal with c = d% = 0.48 and C2 = —0.6 (as predicted in the MCHM5 for £ = 0.3), 
the coupling C2 can be measured, up to a discrete ambiguity, with a precision of ~ 15 — 20% (see 
Fig.§. 

Our partonic analysis of pp — >• hh —> bb^fj should be considered as a first estimate of the LHC 
potentiality, although we expect it to be robust and moderately conservative. For example, we have 
followed a cut-based strategy to reduce the background, although a realistic analysis will certainly 
make use of shape variables and extract the background from data, similarly to what has been done 
by the ATLAS and CMS collaborations in single Higgs searches. Also notice that we did not make 
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use of the information on the total invariant mass distribution, da / dm(bbjj) , which was instead 
used in Ref. [26 1 to further increase the signal significance. Finally, a more precise assessment of 
the LHC sensitivity will require full inclusion of showering and hadronization effects, as well as a 
detector simulation. 
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A Appendix 

We report here the details of the statistical analysis used to derive our results. We follow the 
Bayesian approach N and construct a posterior probability for the total event rate r, 

L(N\r£)ir(r) . 

for given number of observed events iV and luminosity C We denote with ir(r) the prior distribution 
and with L(N\rC) the likelihood function, which we take to be a Poisson distribution 

-rC / C) N 

L{N\rC) = e ^> . (17) 

For the plots of Figs. [§ [7] and [8] we use a prior distribution which is flat for r > and vanishing 
otherwise. The discovery contours of Figs. [6] and [7] are obtained by setting the number of observed 
events to the total value N = (r s + r^)jC expected in each point of the Higgs coupling's parameter 
space: r s = r s (a, c, C2, cfo) for Fig. [6| r s = r s (£) for Fig. [7| We define a point in this space to be 
'discoverable' at a certain luminosity C if the probability of having a total number of events smaller 
than or equal to r^C is below 1%, 

"dr' Pc (r'\(r s +r b )C) < 0.01, (18) 



7 See for example Ref. 45 for a primer. 
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and if the number of observed events (r s + r^)C is 5 or larger. The discovery luminosity is thus 
defined to be the smallest value of C which satisfies these two conditions. 

The plot of Fig. [9] has been derived by expressing the posterior probability as a function of £, 
(pc(£\N) = Pc( r s(£)~\~ r b\N)), and choosing a flat prior for < £ < 1 (vanishing otherwise). For each 
value £th the number of observed events has been set to the total expected value N = (r s (£t/ l )+r&)£. 
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